NASA Technical Memorandum 103726 
AIAA-90-5240 


I-W4I 

£5l<°o 


Experimental Investigation of a Single Flush- 
Mounted Hypermixing Nozzle 


David O. Davis, Warren R. Hingst and A. Robert Porro 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

Second International Aerospace Planes Conference 

sponsored by the the American Institute of Aeronautics and Astronautics 

Orlando, Florida, October 29-31, 1990 


NASA 


EXPERIMENTAL INVESTIGATION OF A SINGLE FLUSH-MOUNTED 

HYPERMIXING NOZZLE 


D.O. Davis* *, W.R. Hingst* and A.R. Pono* 
NASA Lewis Research Center, Cleveland, Ohio 


Abstract 

Reported herein are the results of an experimental 
wind tunnel investigation of a circular supersonic jet (My 
= 3.47) injected at a 10 degree angle into a supersonic 
freestream. The jet penetrates a boundary layer, which has 
a thickness approximately the same as the jet nozzle exit 
diameter. Measurements were made for nominal freestream 
Mach numbers of 1.6, 2.0, 2.5 and 3.0. Three jet total pres- 
sures were run at each freestream Mach number, resulting 
in twelve separate operating conditions. Mean data accu- 
mulated by means of static and total pressure probe in- 
strumentation Eire presented at two axial stations: seven 
jet nozzle diameters upstream and 15 jet nozzle diame- 
ters downstream from where the centerline of the nozzle 
intersects the wind tunnel wall. For one condition at each 
freestream Mach number, the jet air was seeded with a hy- 
drocarbon trace gas and the flow was sampled at the down- 
stream measurement plane to quantify the mean mixing of 
the two streams. Surface oil flow visualization was also used 
to investigate the flow interaction. All results are for air- 
to-air mixing. The measurements indicate the presence of 
two pairs of contra-rotating vortices. One pair follows the 
jet trajectory and tends to split the jet into two streams. 
A smaller pair, rotating in an opposite sense, develops in 
the near wall region. Reported results include Mach num- 
ber and volume fraction distributions in the cross plane, as 
well as jet penetration and mixing efficiency. 


Nomenclature 



= aurea of mixed gas at measurement plane 

Atotal 

= area of wind tunnel at measurement plane 

D 

= nozzle exit diameter (24.2 mm) 

h 

= jet penetration height 

h 12 

= boundau-y layer shape factor (£ 1 / 62 ) 

M 

= Mach number 

N 

= number of moles 

M 

= Mach number ratio ( Mj/Moo ) 

P 

= static pressure 

Pe 6 

= effective back pressure 

Pt 

= total pressure 

Re 

= unit Reynolds number (U 0 0 /^) 
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T t = total temperature 

q = dynamic pressure ratio ( pU 2 )j/(pU 2 )oo 

u' = RMS ELxiad fluctuating velocity component 

U = axial velocity component 

v = volume (mole) fraction 

x, y, z = cairtesian coordinate system (Fig. 2) 

6 = boundary layer thickness 

<$! = displacement thickness 

62 = momentum thickness 

6 = injection angle ( 10 °) 

A = mass flux ratio ( pU)j/(pU)oo 

\i — molecular viscosity 

v — kinematic viscosity (/i/p) 

p = density 

Subscripts 

j = jet condition 

00 = freestream (tunnel) condition 

e = boundary layer edge condition 

Introduction 

An experimental study is being conducted at NASA 
Lewis Research Center to determine the effectiveness of 
various hypermixing nozzle concepts for NASP propulsion 
applications. A single flush nozzle has been selected as 
a baseline model to investigate jet penetration and mixing 
over a rarnge of operating conditions. Results obtained with 
this nozzle will provide a comparison for future testing of 
more 'complicated multi-nozzle configurations (see Fig. 1). 
These results will also provide a data base for Computa- 
tional Fluid Dynamics (CFD) code development. 

Nozzles with shallow injection angles represent a com- 
promise between the extremes of normal injection amd tan- 
gential slot injecton. Jets injected normal to a supersonic 
stream provide good penetration and mixing, but at the 
expense of large total pressure loss. Tangential slot injec- 
tion, on the other hand, contributes all of it’s momentum 
to axial thrust, but with no penetration and little mixing. 
Multiple shallow angle injectors which are ramped and sep- 
arated by expansion slots utilize oblique shock interaction 
and streaunwise vorticity generation to enhance mixing . 1 A 
typical injector block of this type is shown in Fig. 1. In or- 
der to evaluate the effectiveness of an injector configuration 
such as this, it is useful to first consider the case of a flush 
mounted nozzle to determine the effect of jet properties 
alone on penetration and mixing. 

The injection of an underexpanded secondary jet at a 
non-zero injection angle into a primary supersonic stream is 
assumed to be a two-stage process . 2,3 The first stage occurs 
in the region of initial penetration, where the jet is acceler- 


1 



Fig. 1. Swept ramp hypermixing nozzle. 


ated and turned in the direction of the primary stream. Lit- 
tle mixing occurs in this region. The second stage is nearly 
coaxial turbulent mixing of the two streams. The presence 
of the jet acts as an obstacle to the primary stream and 
generates an interaction shock that increases the effective 
back pressure of the jet. The flow expands in the presence 
of the effective back pressure until a normal shock (Mach 
disk) occurs. Schetz et al. A proposed that an analogy ex- 
ists between injection into a supersonic cross stream and 
injection into a quiescent medium. By suitable definition 
of the effective back pressure, the penetration of the Mach 
disk could be predicted from the well-understood theory 
of the latter case. Schetz found that for normal injection 
of a sonic jet 5 and supersonic jet 6 into a supersonic cross 
stream, the penetration of the Mach disk correlated well 
when the effective back pressure was defined as 80% of the 
static pressure behind a normal shock in the freestream. 
Recently, Mays et al. 7 suggested that for shallow angle in- 
jection, the effective back pressure could be approximated 
by assuming that the jet appeared to the freestream as a 
half cone with a half-angle equal to the injection angle. 

Analytic jet penetration models for normal injec- 
tion 3,5,8 and for angled injection 9-11 have been developed. 
These models, however, were formulated to predict the pen- 
etration of the Mach disk and axe restricted to the nearfield 
interaction region. In the downstream region, the jet de- 
cays by turbulent mixing and is distorted by streamwise 
vortices as illustrated in Fig. 2. 12 Correlations for penetra- 
tion in the downstream region have been proposed by Orth 
et al. 2 and Faucher et a/. 13 for normal sonic injection and 
by McClinton 14 for angled sonic injection of multiple jets. 
Rogers 15,16 attempted to correlate his normal single and 
multiple sonic jet data with the model given in Refs. 2 and 



Fig. 2. Sketch of normal jet flowfleld. 12 


13 with limited success. Mays et al 7 found poor agree- 
ment between their shallow angle sonic injection and the 
correlation proposed by McClinton. 

At present, jet penetration and mixing correlations for 
shallow angle supersonic injection into a supersonic stream 
axe virtually non-existent due to the lack of experimental 
data. In addition, detailed data in the transverse plane is 
required to support CFD code validation and calibration. 
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The present study was initiated to help fill this void. The 
nozzle configuration chosen for this study consists of a cir- 
cular supersonic jet (My =3.47) injected at a 10 degree angle 
into a supersonic freestream. The jet penetrates a bound- 
ary layer whose thickness is approximately the same as the 
nozzle exit diameter. 

Experimental Program 

Flow Facility 

The flow facility used in this study was the NASA- 
Lewis Research Center 1 x 1-foot Supersonic Wind Tunnel. 
This wind tunnel is a continuous, open-loop flow facility 
with a Mach number range from 1.3 to 5.5 provided by in- 
terchangeable nozzle blocks. A schematic of the wind tun- 
nel test section with the injector nozzle installed is shown in 
Fig. 3. The nozzle is located halfway across one of the wind 
tunnel sidewalls whose height is 30.5 cm. A splitter plate, 
which supports the probe traverse mechanism, is located at 
the midspan of the wind tunnel. Reference coordinates for 
this configuration and the two axial locations where data 
were accumulated are also indicated in Fig. 3. 

Pressure Probe Instrumentation 

The test section was outfitted with both static and 
total pressure instrumentation. Static pressure at the up- 
stream station was sampled by means of 0.51 mm diameter 
wall taps. Boundary layer surveys at the upstream location 
were measured with a flattened pitot tube having outer di- 
mensions of 0.81 x 0.41 mm at the tip. Static pressure in 
the cross plane at the downstream station was measured 
with a 15 tube cone-cylinder probe rake. The half-angle of 
the cone was 5°, the diameter of each probe was 1.02 mm 
and the tubes were spaced 9.53 mm apart. Pitot pressure 
at the downstream station was measured with a 15 tube 
pitot rake, which had the same nominal dimensions as the 
static pressure rake. The rakes were positioned such that 
the probes lay in a plane parallel to the i-y plane. 


Trace Gas Instrumentation 

A hydrocarbon trace gas technique was used to mea- 
sure the mean mixing of the two streams. The supply line 
to the nozzle was seeded with ethylene at a location ap- 
proximately 60 supply line pipe diameters (Pipe I.D.=5.08 
cm) upstream from' the injector plenum to ensure thor- 
ough mixing. An Edwards Model 825 mass flow controller 
maintained the ethylene flow at a constant rate. At the 
downstream measurement plane, the flow was sampled by 
means of the pitot tube rake described in the previous sec- 
tion. The sampled gas passed through a series of pumps 
and back pressure regulators before entering a Gow-Mac 
Model 23-500 totad hydrocarbon analyzer. Because the gas 
route must be completely purged of the previous data point 
sample before a current reading is taken, the time between 
readings is relatively long; approximately 30 seconds. In 
contrast, a pressure measurement typically taices 5-10 sec- 
onds. The output of the analyzer is a voltage linearly pro- 
portional to the number of carbon atoms flowing through 
the unit. The anailyzer is calibrated with zero gas adr and 
an NBS grade reference mixture of propame amd zero gas air 
(475 ppm). For the two jet total pressures considered for 
the trace gas measurements, namely; 2.07 amd 2.76 MPa, 
the ethylene concentration in the jet plenum was observed 
to be 150 and 112 ppm, respectively. 

When properly calibrated, the traice gas analyzer dis- 
plays the volume (mole) fraction of ethylene, which for the 
present case is equivalent to the mass fraction. Readings 
taLken in the measurement plane, when normadized by the 
volume fraction in the jet plenum, yield the volume fraction 
of the jet to tunnel flow: 

(-^ tthyltnt / N total) Nj / ,\ 

Vi ~ (N.tk,Un ,/Nj) ~ Ntou, 1 1 
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Table 1. Operating Conditions 


Case 

Moo 

P,,oo (kPa) 

Ptj/Pt.OO P j / P oo 

Pi/p* 

Uj/U oo 

A 

7 

1 

1.6 

89.6 

15.4 

0.86 

0.69 

1.45 

2.83 

4.10 

2 

1.6 

89.6 

23.1 

1.29 

1.03 

1.45 

4.25 

6.15 

3* 

1.6 

89.6 

30.8 

1.72 

1.37 

1.45 

5.67 

8.20 

4 

2.0 

137.8 

10.0 

1.03 

0.77 

1.26 

2.48 

3.14 

5* 

2.0 

137.8 

15.0 

1.54 

1.15 

1.26 

3.73 

4.71 

6 

2.0 

137.8 

20.0 

2.05 

1.53 

1.26 

4.97 

6.28 

7 

2.5 

172.4 

8.00 

1.79 

1.21 

1.13 

3.11 

3.51 

8* 

2.5 

172.4 

12.0 

2.69 

1.81 

1.13 

4.66 

5.27 

9 

2.5 

172.4 

16.0 

3.58 

2.42 

1.13 

6.21 

7.03 

10 

3.0 

206.8 

6.67 

3.21 

1.95 

1.05 

4.16 

4.37 

11* 

3.0 

206.8 

10.0 

4.82 

2.93 

1.05 

6.24 

6.55 

12 

3.0 

206.8 

13.3 

6.42 

3.90 

1.05 

8.32 

8.74 

M,= 

3.47, T tJ 

=T t>oo =294.0 K, 0=10° 

(all cases) 






Results and Discussion 
Operating Conditions 

Data were accumulated for frees tram Mach numbers of 
1.6, 2.0, 2.5 and 3.0. At each Mach number, three jet total 
pressures were considered, resulting in a total of 12 operat- 
ing conditions. However, due to the long time required to 
analyze a trace gas sample, trace gas measurements were 
made at only one condition per Mach number. The 12 op- 
erating conditions are summarized in Table 1. Asterisks by 
the Case number indicate that trace gas measurements were 
made. The total temperature of the streams wets nominally 
the same (T*=294° K). 

Following Mays et al? , the effective back pressure (P c b) 
is defined as the surface pressure on a cone in the freestream 
which has a half- angle equal to the injection angle. Assum- 
ing this to be a valid approximation, then test Cases 1 and 
4 axe at a overexpanded condition. The pressure ratio for 
Case 2 is nearly at a matched condition, while the remain- 
der of the cases are at underexpanded conditions. 

Disturbances in the freestream generated by the pres- 
ence of the jet will propagate outward at an angle slightly 
larger than the Mach angle associated with the freestream 
condition. These waves will reflect from the wind tunnel 
walls and pass through the jet stream. This is an im- 
portant consideration because the interaction between an 
oblique shock and a mixing layer is known to increase the 
rate of mixing between the streams 1,17 . The path of these 
disturbances was estimated for the present configuration. 
For the two lower freestream Mach numbers, multiple re- 
flections occur before crossing the measurement plane. For 
the two higher Mach numbers, the waves pass through the 
measurement plane without crossing the jet stream. 

Boundary Laver Profiles 

Inasmuch as this study is intended to support concur- 
rent CFD efforts at the NASA Lewis Research Center, it 
was considered important to provide well documented up- 
stream flow conditions. To this end, boundary layer profiles 
were measured along the wind tunnel wadi bisector at a loca- 
tion seven nozzle diameters upstream from where the nozzle 
centerline intersects the wadi (x=0). Meam velocity profiles 


for the four freestream operating conditions are plotted in 
Fig. 4. Previous hot-wire measurements and anadysis of 
the profiles have shown that the boundary layers axe fully 
turbulent. Integral parameters as well as freestream turbu- 
lence intensity axe summarized in Table. 2. 


Table 2. Upstream Flow Condition 


Moo 

1.6 

2.0 

2.5 

3.0 

6/D 

0.924 

1.081 

1.265 

1.110 

6,/D 

0.175 

0.230 

0.290 

0.347 

St/D 

'0.078 

0.083 

0.081 

0.075 

H n 

2.240 

2.770 

3.580 

4.630 

U'/U oo (%) 

0.470 

0.450 

0.380 

0.410 


D = 24.2 mm 



Fig. 4. Mean velocity profiles at x/D=- 7. 
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Flow Visualization 

A qualitative assesment of the near wall flow behav- 
ior was obtained by means of surface oil flow visualization. 
The entire wall of the test section was painted with a mix- 
ture of SAE HOW oil and a flourescent dye. The oil was 
brushed on with strokes perpendicular to the primary flow 
direction. The wind tunnel was run until a steady pat- 
tern was established (approximately 15 minutes) and then 
rapidly shutdown to preserve the flow pattern. An ultra- 
violet light was used to illuminate the dye pattern. The 
results for Case 8 axe shown in Fig. 5. The results for the 
other flow conditions exhibited the same general behavior 
as this case. The flow pattern shows no indication of flow 
separation, although a small region may exist immediately 
downstream of the nozzle. In general, the flow is displaced 
around the jet and then converges back towards the cen- 
terline. The structure downstream of the nozzle is a result 
of a pair of contra-rotating vortices that have the common 
flow towards the wall (-y direction). Oil in this region is 
displaced away from the centerline of the wind tunnel. Nu- 
merical predictions performed at the NASA Lewis Research 
Center using the RPLUS 18 flow solver indicate the presence 
of these vortices. This feature of shallow angle injection is 
believed to be previously unreported. 



Mach Contours 

The static and pitot pressure surveys were used to cal- 
culate Mach number distributions in the cross plane. These 
results are shown for the twelve operating conditions in 
Figs. 6a through 61, respectively. The peak Mach number 
observed in the cross section is also shown in the figures. A 
physical limitation of the probe traversing mechanism only 
allowed the surveys to be extended to the z/D =- 2 loca- 
tion. This limitation was corrected for the trace gas mea- 
surements. In general, the flowfield for all cases is charac- 
terized by the pressure induced vortex pair associated with 
a jet in a cross flow. 12 These vortices create the outwaxd 
bulging of the contours away from the wall and axe in addi- 
tion to the aforementioned neax wall vortex pair that were 
evident in the oil flow visualization (Fig. 5). Although the 
near wall vortices are not resolved by the transverse flow- 
field measurements, evidence of their presence can be seen 
in Figs. 6i and 61, where the neax wall contour (M=1.4) in 
the vicinity of z/D=0 bulge inward towards the wall. A 
certain degree of asymmetry is observed to be present at 
the lower freestream Mach numbers (Figs. 6a through 6f). 
A similar asymmetry was observed for the 15° sonic helium 
jet investigated by Mays ti a/. 7 , which they attributed to 
the high streamwise vorticity in the flow. For Cases 8-12, 
the effective pressure ratio causes sufficient expansion of the 
jet flow, producing local regions in the measurement cross 
plane where the Mach number exceeds the nozzle exit Mach 
number. It should be noted here that for Mach numbers 
greater than approximately 4.5, the operating total tem- 
perature of the streams ( T t = 294° K) was not high enough 
to prevent condensation of the air stream. The number of 
Mach number peaks in the cross plane tends to increase as 
the freestream Mach number increases. At Moo=1.6 and 2.0 
(Figs. 6a through 6f), two peaks in the Mach number are 
observed which are nominally symmetric about the y-axis. 
For Moo=2.5 (Figs. 6g through 6i), a third peak appeaxs 
along the plane of symmetry (y-axis). For Moo=3.0 (Figs. 
6j through 61), in addition to the peak along the y-axis, two 
pair symmetrically located about the y-axis axe observed. 


Fig. 5. Oil flow visualization for Case 8. 
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z/D z/D 


c) Case 3, Moo = 1.6, q=8. 20. 


f) Case 6, Moo=2.0, £=6.28. 


Fig. 6. Mach number contours at x/D=15. 
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z/D z/D 


g) Case 7, M co = 2.5, g=3.51. 


j) Case 10, Moo=3.0, g=4.37. 




z/D z/D 

h) Case 8, M a 0 =2.5, q= 5.27. k) Case 11, M oo =3.0, g=6.55. 




2 /D z/D 

i) Case 9, A/*, =2.5, $=7.03. 1) Case 12, M oo =3.0, «=8.74. 


Fig. 6. Mach number contours at x/D= 15. 
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Mixing Behavior 

Trace gas measurements were made for Cases 3, 5, 
8 and 11. Contours of the volume (mole) fraction in the 
measurement plane, as well as the peak concentration mea- 
sured, are shown in Figs. 7a through 7d. The vortex in- 
duced pattern of the contours is similar to that observed in 
the Mach number contours. The asymmetry observed for 
the lower freestream Mach numbers is also present in the 
trace gas results. The spreading of the jet is observed to be 
largest for Case 3. A feature of this case that is not observed 
in the others is the outward bulging of the lowest contour 
level in the vicinity of y/D= 1, z/D=± 2. This is believed to 
be a result of the interaction between the reflected oblique 
shock sind the jet stream. It is also likely that this inter- 
action enhances mixing for this case. It is interesting to 
note that, with the exception of Case 3, the spreading of 


the jet increases as the freestream Mach number increases 
( Mj/Moo decreases), which is opposite to what would be 
expected with ail other variables being constant. 

In order to quantify the mixing between the two 
streams, the mixing efficiency proposed by Hollo ti a/. 17 
will be used. The mixing efficiency is defined as the per- 
centage of the test section area that is considered to be 
“mixed”. Hollo suggests that the flow is considered mixed 
when the volume fraction (v) is between the static flamma- 
bility limits for hydrogen in air, namely; 4 to 74%. The 
choice of this range is somewhat arbitrary, but is useful 
for comparative purposes. For the present study, the lower 
limit was specified as the lower flammability limit for hy- 
drogen in air (t>£,=0.04), but in order to get a better idea 
of the distribution of the mixture, mixing efficiencies based 
several upper limits are presented (vc/=0.25, 0.5, 0.74 and 




z/D z/D 


a) Case 3, Moo = 1.6, <?=8. 20. 


c) Case 8, Moo=2.5, g=5.27. 




z/D z/D 

b) Case 5, Moo=2.0, £=4.71. d) Case 11, Moo=3.0, £=6.55. 


Fig. 7. Volume fraction contours at z/D= 15. 
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1.0). The mixing efficiency defined in this manner is plot- 
ted against the dynamic pressure ratio in Fig. 8. The curve 
labeled V(/ = 1.0 represents the fraction of the test section 
area through which at least a trace (t> > 0.004) of the jet 
fluid passes and is an indication of the degree of spread- 
ing of the jet. The curve labeled vy= 0.74 represents the 
fraction of the test section area that is within the static 
flammability limits for hydrogen in air. For Cases 3 and 5, 
essentially all of the mixture in the cross plane is below the 
upper flammability limit (see Figs. 7a and 7b). Although 
the data is limited, the spreading of the jet appears to cor- 
relate strongly with the dynamic pressure ratio. Another 
interesting feature is that the cases which had the lowest 
(Case 5) and highest (Case 11) peak concentration levels, 
also had the smallest and second largest mixing efficiencies, 
respectively. This implies that peak concentration is not a 
good indicator of overall mixing. 



Fig. 8. Mixing efficiency vs. q at x/D= 15. 


Jet Penetration 

For the aerodynamic measurements, the jet penetra- 
tion is defined as the location along the y-axis (symmetry 
plane) where the Mach number is at a peak. Similarly, for 
the trace gas measurements, the penetration is defined as 
the location along the plane of symmetry where the peak 
concentration occurs. The penetration heights as defined 
above are plotted versus dynamic pressure ratio in Fig. 9. 
At lower values of y, the location of peak concentration cor- 
responds with the peak Mach number location, but then 
deviates at higher values of q. The data, however, are too 
limited to draw any definitive conclusions. For the three 
lower freestream Mach numbers, the slope of the penetra- 
tion curves are seen to be similar, but a decrease in sensitiv- 
ity to dynamic pressure ratio is observed for the Afoo =3.0 
data. 



q 

Fig. 9. Jet penetration vs. q at x/D=15. 


Penetration correlations for normal sonic injection 
have been proposed and sure reviewed by Rogers. 15 These 
correlatons are typically of the form: 

h/D = a(q) b (x/D) c (2) 

The present data indicate a Mach number ratio dependence 
so that an equation of the form: 

h/D = a{q) i {M) c {x/ D) 11 (3) 

was used to attempt to correlate the aerodynamic penetra- 
tion data with dynamic pressure ratio and Mach number 
ratio. A least-squares curve fit of the data to this equation 
is shown in Fig. 10. At least over the conditions considered 
in the present study, this relation provides a reasonable de- 
scription of the jet penetration. 



Fig. 10. Jet penetration correlation. 
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Concluding Remarks 

Flowfield measurements in a plane 15 nozzle diame- 
ters downstream of a supersonic jet injected at a 10° angle 
into a supersonic stream have been presented. The experi- 
mental results, along with the supporting CFD calculations, 
show that the flowfield is characterized by two vortex pairs. 
The pair furthest from the wall significantly distorts the jet 
mean flowfield. The results also indicate that jet penetra- 
tion and mixing correlate strongly with dynamic pressure 
ratio. In addition, at a given dynamic pressure ratio, the 
jet penetration was observed to be inversely proportional to 
Mach number ratio. A correlation has been presented that 
describes the jet penetration for the present configuration, 
but additional data at other axial locations is needed for 
verification. 

The present data set can be used for CFD calibration 
and validation purposes, inasmuch as cross plane data is 
presented over a fairly wide range of operating conditions 
and upstream conditions are well defined. For the jet flow, 
calculations should begin at the throat of the nozzle since 
it is unlikely that the flow is uniform at the jet exit. The 
splitter plate, as well as the wind tunnel walls, should be 
modeled to account for oblique shock wave interaction in 
the flowfield. 
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